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ABSTRACT 


In-bore  spectroscopic  observations  of  a  free-flowing  plasma  for 
wavelengths  between  3000°A  and  6250°A  were  made  for  a  series  of  railgun 
firings.  In  this  report  the  experimental  arrangement  for  the  series  is 
described  dand  the  results  are  presented.  Atomic  species  present  In  the 
plasma  are  identified  from  absorption  and  emission  spectra.  Calculations  of 
the  degree  of  Ionisation  of  these  species  are  then  used  to  produce  a 
temperature  estimate  between  11  x  103  K  to  25  x  103  k  for  the  railgun 
plasma.  A  comparison  of  theoretical  values  of  the  plasma  electrical 
conductivity  with  values  obtained  from  muzzle  voltage  records  Indicates  that 
the  peak  plasma  temperature  Is  approximately  11  x  103  K,  which  Is  consistent 
with  estimates  obtained  from  the  spectra. 
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cross-sectional  area  of  railgun  bore,  m 

ratio  of  the  number  density  of  electrons  to  the  number  density  of 
heavy  particles 

magnetic  field,  T 

ratio  of  conductivity  to  that  in  a  Lorentz  gas 

degeneracy  factor  of  an  ion  in  the  i-th  excited  level  of  the  j-th 
ionisation  state 

current  flowing  in  railgun,  A 

ionisation  energy  needed  to  ionise  the  atom  j  times 
current  density,  A/m2 

dimensionless  parameter  in  Saha's  equation 
Boltzmann's  constant,  1.3807  x  io-23  J/K 
wavelength  of  light.  Angstroms 
inductance  per  unit  length  of  railgun,  H/m 
electron  mass,  kg 

atomic  mass  of  ionised  or  neutral  atom,  kg 

electron  number  density 

dimensionless  distance 

plasma  pressure,  Pa 

plasma  mass  density,  kg/m3 

plasma  temperature,  K 

energy  of  the  1-th  electronic  level  of  an  ion  in  the  j-th  ionisation 
state 

ratio  of  the  number  density  of  ions  ionised  j  times  to  the  total 
number  density  of  heavy  particles 

electronic  partition  function 


ionic  charge 
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TEMPERATURE  ESTIMATES  FOR  A  FREE 
FLOWING  RAILGUN  PLASMA 


IXTRODUCTIOS 


Thi3  report  documents  spectroscopic  measurements  for  wavelengths 
between  3000°A  and  6260°A  of  a  free  flowing  plasma  taken  In  a  series  of 
firings  using  a  sma  1  bore  ( 10  mm  x  10  mm)  rallgun  and  reports  estimates  of 
plasma  temperature  derived  from  these  observations.  The  observations  were 
made  along  the  bore  of  the  rallgun  using  a  0.75  m  Jarrel  Ash  Spectrometer  and 
recorded  on  photographic  film.  The  spectrum  is  typically  a  continuum  with  a 
number  of  absorption  lines. 


The  plasma  armatures  were  created  by  electrically  exploding  a  number 
of  different  metallic  foils.  Various  Ionised  elements  within  the  plasma  were 
Identified  from  the  spectra.  Calculations  for  the  degree  of  ionisation  for 
first-  and  second-ionised  species  were  then  used  to  obtain  an  estimate  of  the 
peak  plasma  temperature.  These  calculations  are  based  upon  certain 
assumptions  which  are  discussed  In  this  report. 

This  report  is  organised  into  7  major  sections,  in  section  1,  a 
brief  history  of  the  experimental  and  theoretical  work  on  estimation  of  the 
temperature  of  rallgun  plasmas  is  presented  whilst  in  Section  2,  the 
experimental  arrangement  for  the  series  Is  described.  In  section  3,  the 
theory  used  to  obtain  the  temperature  estimate  from  the  spectra  is  presented 
in  detail.  The  experimental  results  are  presented  and  discussed  in 
Section  4.  The  results  of  section  4  are  compared  with  the  results  obtained  by 
other  rallgun  workers  In  Section  5.  In  Section  6,  the  temperature  is 
estimated  from  the  electrical  conductivity  of  the  plasma  and  compared  with  the 
experimental  estimates  of  the  temperature  reported  in  Section  4.  Finally,  in 
section  7  we  present  a  conclusion  from  the  results  described  in  this  report. 
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1. 


BACKGROUND 


In  1978,  Rashleigh  and  Marshall  III  reported  details  of  experiments 
with  an  inductively-driven  rallgun  In  which  a  cubic  lexan  projectile  was 
accelerated  to  5.9  Ian  s-1.  This  pioneering  work  has  been  followed  by 
considerable  rallgun  development  and  related  equation  of  state  research  on 
plasmas.  However,  little  experimental  research  into  the  tempeiature  of  the 
plasma  armature  within  the  bore  of  the  rallgun  has  been  done. 


The  first  attempt  to  measure  the  temperature  of  a  rallgun  plasma  In 
recent  times  was  made  by  Parkes  and  Strachan  (21.  Essentially  they  recorded 
the  spectral  emissions  from  the  muzzle  flash  after  the  projectile  had  left  the 
railgun.  By  identifying  ions  in  the  muzzle  flash,  an  estimate  of  the  peak 
temperature  was  made  by  assuming  that  1/5  of  the  Ionisation  energy  of  an  Ion 
In  a  plasma  represented  a  measure  of  its  temperature.  For  their  experiment,  a 
temperature  of  2  58  x  103  K  was  estimated. 


in  1980,  McHab  ( 3 1  calculated  the  temperature  of  the  plasma  in  the 
Rashlelgh-Marshall  rallgun  by  using  a  simple  form  of  Saha's  equation  for  the 
degree  of  ionisation.  Assuming  that  only  single  ionisation  occurred  in  the 
plasma,  he  predicted  that  the  plasma  would  have  a  temperature  of  (44  ±  13)  x 
103  K  and  an  ion  density  of  <8  ±  l)  x  1025  m-3  for  a  current  of  300  kA  in  the 
circuit. 


Powell  and  Batteh  Ml  (hereafter  referred  to  as  PB)  followed  with 
more  rigorous  one-dimensional  magnetohydrodynamlc  (mhd)  calculations  and  a 
form  of  saha's  equation  which  allowed  for  the  occurrence  of  double  ionisation 
in  the  plasma.  Some  major  assumptions  In  their  model  were  that  the  plasma 
obeyed  the  Ideal  gas  law  and  that  the  plasma  was  in  local  thermodynamic 
equilibrium  ( LTE) •  when  applied  to  the  Rashleigh  and  Marshall  (hereafter 
referred  to  as  RM)  experiment  111,  their  code  estimated  plasma  temperatures 
peaking  at  72  x  103  K  with  electron  densities  of  2.2  x  1026  m~3  at  peak 
pressures  around  230  MPa  and  peak  currents  of  300  kA.  This  work  was  later 
developed  into  a  two-dimensional  ( 2— D>  code  by  Powell  (51,  which  made  similar 
assumptions  to  those  In  the  1-D  model.  Estimates  of  peak  temperatures  in  the 
RM  experiment  were  revised  to  47  x  I03  K  for  electron  densities  of  2.5  x  1026 
m-3.  Powell  accounted  for  the  temperature  drop  as  being  due  to  the  greater 
surface  area  available  for  radiation  of  energy  from  the  plasma. 


In  1983,  Thlo  161  published  details  of  his  1-D  railgun  computer 
simulation  code  in  which  he  assumed  only  single  ionisation.  His  code,  PARA, 
used  quasl-statlc  MHD  equations  to  model  the  time  varying  properties  of  the 
plasma  whereas  McNab(3!  and  PB  (4)  were  essentially  steady  state  analyses. 
Thlo  predicted  temperatures  of  26  x  103  K  at  50  MPa  for  peak  currents  of  125 
kA  in  a  6  x  8  mm  bore  rallgun. 


During  1983,  spectroscopic  observations  o,.  the  plasma  armature  were 
made  on  a  small  bore  rallgun  at  MRL  -  the  RAPID  railgun.  Spectroscopic 
observations  were  made  by  Richardson  and  Clark  (71  of  a  puff  of  plasma  ejected 
through  a  hole  in  the  sidewall  of  a  rallgun.  This  work  revealed  that  the 
plasma  was  optically  dense  because  the  spectrum  contained  numerous  line 
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reversals  superimposed  upon  a  continuum  background.  They  suggested  that  more 
useful  spectroscopic  results  would  be  obtained  In  the  ultra  violet  (irv) 
spectral  region  where  the  continuum  would  reduce  In  intensity. 

The  above  work  was  continued  by  Richardson  In  1983  with 
spectroscopic  studies  of  the  muzzle  flash  of  a  RAPID  rallgun  C a ) .  Thi3  work 
revealed  that  ablation  of  rail  material  into  the  plasma  armature  occurred. 
Using  his  spectral  results,  Richardson  estimated  that  the  plasma  temperature 
was  at  least  30  x  103  K.  Peak  current  in  these  rallgun  experiments  was 
approximately  80  kA. 


Also  In  1983,  Marshall  (91  obtained  an  estimate  of  the  plasma 
armature  temperature  of  43  x  103  K  for  a  plasma  current  of  approximately  30  kA 
by  treating  the  flow  of  a  plasma  through  a  hole  in  the  3ldewall  of  a  rallgun 
as  the  flow  of  a  neutral  gas  through  a  nozzle. 


Gathers  and  Hord  (101  made  measurements  of  rallgun  plasma 
temperatures  by  using  two-colour  pyrometery  and  assuming  black  body  radiation. 
They  reported  a  mean  peak  temperature  of  approximately  5.1  x  lo3  K  for  a 
stationary  arc  in  a  rallgun.  The  magnitude  of  the  current  flowing  through  the 
stationary  arc  was  not  reported. 


In  1986  at  the  Electromagnetic  Launcher  Technology  Conference, 
Marshall  till,  Tidman  et  al  (121  and  Uglum  1131  presented  papers  on  rallgun 
plasmas.  Marshall  (111  reported  a  temperature  of  50  x  103  K  in  a  plasma 
emerging  from  the  muzzle  of  a  rallgun.  The  temperature  was  calculated  from  a 
measurement  of  the  speed  of  sound  In  the  plasma.  Peak  current  in  the  railgun 
was  approximately  loo  kA. 

using  a  turbulent  hydrodynamic  approach,  Tidman  et  al  (12) 
calculated  peak  temperatures  less  than  25  x  103  K  for  the  plasma  surface.  In 
these  calculations  the  current  and  bore  of  the  rallgun  were  respectively  200 
kA  and  lo  mm  in  diameter.  According  to  these  authors,  turbulence  causes  a 
cooling  of  a  rallgun  plasma  with  a  resultant  Increase  In  Its  resistivity. 


Finally,  Uglum  (131  reported  spectroscopic  results  of  a  rallgun 
muzzle  flash.  This  work  was  similar  to  that  done  by  Parkes  et  al  (2). 

Spectral  lines  belonging  to  cul.  Cull,  All  and  sll  were  identified. 

As  Is  evident  from  the  above  history,  little  work  in  the  field  of 
experimental  plasma  diagnostics  for  the  temperature  measurement  of  a  rallgun 
plasma  armature  within  the  bore  of  the  rallgun  has  been  reported.  The  plasma 

armature  Is  not  an  Ideal  plasma  to  model  for  it  is  composed  of  many  different 
ionic  species  (8i  and  on  some  occasions  its  interior  has  density  and 
temperature  gradients  (5,14).  Accordingly,  numerous  assumptions  have  to  be 
made  about  the  plasma  so  that  theory  can  yield  some  insight  Into  its 
properties.  These  are  discussed  later  In  this  report. 


The  methods  for  estimating  the  plasma  temperature  presented  In  this 
report  should  be  applicable  to  all  the  experiments  discussed  In  this  section, 
even  though  rallgun  currents  were  varied. 
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2 .  EXPERIMENTAL  ARRANGEMENT 

The  experimental  arrangement  of  the  MRL  rallgun  la  shown  in 
Figure  1.  A  200  * F,  10  XV  capacitor  (Maxwell,  Model  33778)  waa  attached  to  a 
l  turn,  1  * H  Inductor  manufactured  from  100  x  6 . 3  mm  copper  bar.  The  main 
switch  was  a  simple  copper-foil  tipped-plunger  which  closed  a  gap  between  two 
copper-tungsten  cylindrical  electrodes.  The  self-activating  crowbar-switch 
was  positioned  directly  between  the  terminals  of  the  capacitor:  it3  mode  of 
operation  is  described  in  Reference  tisl.  Each  firing  was  carried  out  with 
the  capacitor  ban*  charged  to  6  XV.  The  peaX  current  was  between  78.1  XA  and 
79.2  xa  for  each  firing. 


figure  1  MRL  rallgun  and  energy  store 


The  rallgun  (which  was  open  to  the  atmosphere)  was  500  mm  long, 
manufactured  from  two  27  mm  thlcX  sheets  of  laminated  polycarbonate  and 
similar  in  design  to  the  earlier  RAPID  rallgun  system  (161.  The  bore  was  10 
mm  square  as  shown  in  Figure  2 .  cadmium-copper  rails  were  used  because  of 
their  high  conductivity  and  tensile  strength.  The  cadmium  content  was  0.6 
percent. 


The  initial  plasma  was  created  from  the  electrical  explosion  of  a 
foil  which  shorted  the  rails  at  the  breech  end  of  the  rallgun.  The  foil  was 
held  in  place  by  wrapping  It  around  a  polycarbonate  plug  which  was  pushed 
between  the  rails  (Figure  3).  A  wedge  behind  the  plug  secured  the  plug  when 
the  foil  was  exploded. 


Closure  of  the  plunger-switch  In  the  circuit  (Figure  4)  produced  an 
electrical  explosion  of  the  metallic  foil,  which  generated  a  plasma 
armature.  The  plasma  was  accelerated  along  the  rails  by  the  Lorsntz  (J  x  B) 
force.  It  should  be  noted  that  there  was  no  projectile  ahead  of  the  plasma  in 
this  series  of  firings  and  hence  the  term  "free  flowing"  is  used. 
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FIGURE  2  Cross-section  of  the  rallgun 


FIGURE  3  Means  of  holding  the  foil  between  rails 


FIGURE  4  Circuit  diagram  of  the  railgun 


The  muzzle  voltage  and  current  were  recorded  by  the  transient 
recorders  described  in  Reference  (16! .  After  each  firing  the  rails  were 
cleaned  by  swabbing  the  bore  with  alcohol. 

Creation  and  acceieratlon  of  the  plasma  armature  along  the  railgun 
bore  was  recorded  photographically  through  the  transparent  wall  of  the  railgun 
by  an  STL  image-converter  camera.  The  camera  was  positioned  at  the  side  of 
the  launcher  so  that  it  viewed  the  entire  bore  length,  (Figure  5).  Marker 
bars  along  the  length  of  the  gun  provided  position  co-ordinates  on  the 
photographic  record.  The  camera  was  used  in  a  streak  mode  with  a  total  streak 
time  of  200  *s  or  400  * s.  Triggering  of  the  camera  was  achieved  by  directing 
the  light  pipe  from  the  optical  trigger  input  at  the  main  arc  switch. 

Light 


FIGURE  5  Relative  positions  of  the  image-converter  camera  and 
the  Jarrell  Ash  spectroscope 


The  Jarrel  Ash  spectroscope  was  positioned  approximately  3  m  away 
and  directly  in  front  of  the  railgun  muzzle  as  shown  in  Figure  S.  The 
entrance  slit  and  optical  axis  of  the  spectroscope  were  aligned  with  the 
lengthwise  central  axis  of  the  railgun-bore.  Fused-silica  (Herasll)  lenses, 
which  are  transparent  to  ultra-violet  (uv)  light,  provided  Illumination  of  the 
entrance  slit  of  the  spectroscope  with  light  emanating  from  the  plasma.  The 
slit  width  was  set  at  20  #m,  5  mm  high,  and  well-focussed  lines  were  recorded 
on  150  mm  x  loo  mm  Ilford  HP5  film.  Because  there  was  low  light  illumination 
at  the  focal  plane,  the  film  had  to  be  force-developed  in  microphen  to  6400 
asa  to  obtain  a  good  photograph  of  the  spectrum. 

Reference  lines  were  obtained  on  the  film  by  using  an  atomic 
absorption,  hollow  cathode,  copper  reference-lamp  placed  along  the  optical 
axis  near  the  muzzle  of  the  railgun.  Adjustable  optical  masks  at  the  focal 
plane  of  the  spectroscope  enabled  the  reference-spectrum  to  overlap  the 
plasma-spectrum  without  movement  of  the  film.  Thus  the  spatial  relationship 
between  the  reference  and  the  plasma  spectrum  was  direct  and  fixed. 


The  recorded  spectrum  was  integrated  over  the  time  interval 
beginning  with  the  creation  of  the  plasma  and  ending  when  there  was  no  current 
In  the  railgun.  The  spectral  range  was  from  3000°A  to  6260°a. 


spectral  lines  were  resolved  on  the  film  record  by  using  a  model  24- 
201  Jarrel  Ash  microphotometer.  Chart  records  of  the  copper  reference- 
spectrum  were  made  from  each  record  in  order  to  determine  the  dispersion 
factor,  (wavelength/unit  distance).  The  copper  lines  chosen  for  this 
determination  were  3247. 54°A  and  5218. 20°A.  The  resultant  wavelength 
determinations  were  accurate  to  better  than  0.5°A  on  the  reference  films.  It 
should  be  noted  that  a  dispersion  factor  was  required  for  each  film  because 
each  film  was  slightly  different.  The  reason  for  these  differences  i3  unknown 
but  is  probably  due  to  shrinkage  or  expansion  of  the  film  during  processing. 


A  different  type  of  metallic  foil  was  used  in  each  firing 
(Table  l).  The  foil  materials,  together  with  the  materials  in  the  launcher 
and  the  atmosphere,  provided  a  range  of  1st  and  2nd  ionisation  potentials  from 
5.21  ev  to  35.12  ev,  (Table  2).  The  choice  of  foils  was  determined  by  the 
requirement  for  a  spread  of  ionisation  energies  as  explained  in  the  next 
section . 


The  masses  of  aluminium  and  zinc  were  chosen  to  provide 
(approximately)  the  same  number  of  atoms  of  each  material.  This  was  not  done 
for  Baci2  because  its  powdered  form  required  wrapping  in  aluminium  foil  when 
packed  between  the  rails. 
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TABLE  1 


Foils  used  in  the  series  of  firings 


Shot  No. 

Foil  Material 

Mass 

1 

Aluminium  Foil 
(o.oso  mm  thick! 

9.3  mg 

2 

Zinc  Foil 
(0.025  mm  thick) 

22.3  mg 

3 

BaCl2.2H20 
(wrapped  in  A1 
foil) 

59  mg  BaCl2.2H2o 

29  mg  A1 

TABLE  2 


List  of  elements  and  their  respective  ionisation  energies  In) 


Neutral  atom  i  ion 

Ionisation  Energy  (ev) 

source 

cu 

7.724 

Cued  rails 

Cu+ 

20 . 29 

CuCd  rails 

Cd 

8 . 991 

Cued  rails 

Cd+ 

16. 904 

CuCd  rails 

A1 

5 .984 

Foil  A1 

Al  + 

18 . 823 

Foil  A1 

Zn 

9 .391 

Foil  Zn 

Zn+ 

17.96 

Foil  Zn 

Ba 

Ba+ 

5.21 

10. 001 

Foil  (Baci2.2H2o  salt) 

Foil  <BaCl2-2H2o  salt) 

Cl 

13,01 

Foil  <BaCl2.2H20  salt) 

Cl  + 

23 . 8 

Foil  (Baci2.2H2o  salt) 

H 

13.595 

Foil  <BaCl2-2H20  salt)  ( 

H  + 

- 

railgun  body 

0 

13.614 

Atmosphere  6  railgun  body 

0  + 

35.108 

Atmosphere  s  railgun  body 

c 

11.26 

Atmosphere  i  railgun  body 

c+ 

24  .  376 

Atmosphere  S  railgun  body 

N 

14.53 

Atmosphere 

N+ 

29 . 593 

Atmosphere 

s 
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THEORY 


From  the  work  of  Richardson  [71,  It  Is  known  that  the  spectrum 
(within  the  visible  region)  produced  by  a  rallgun  plasma  consists  of  line 
reversals  superimposed  upon  a  continuum.  Line  reversals  indicate  that  the 
plasma  is  optically  thick  (strongly  absorbing). 

Cowan  and  Dieke  118)  have  shown  that  line  reversal  occurs  when  the 
ratio  of  the  number  density  of  emitting  atoms  to  the  number  density  of 
absorbing  atoms  is  not  constant.  This  occurs  when  a  source  is  not  uniformly 
excited.  For  such  sources,  eg.  metal  arcs  and  most  likely  railgun  plasmas,  it 
is  usual  for  the  excitation  to  decrease  from  some  inner  portion  of  the  source 
to  the  edge  of  the  source,  [181.  Furthermore,  if  there  is  only  a  thin,  cool 
(absorbing)  layer  at  the  edge  of  the  source  and  excitation  increases  rapidly 
towards  the  interior,  it  is  highly  possible  that  in  such  cases,  line  reversals 
can  be  used  to  identify  the  excited  atoms  in  the  interior  of  the  plasma,  in 
addition,  according  to  Reference  (191,  line  reversals  will  occur  for  ionised 
species.  It  is  assumed  that  the  above  argument  for  excited  atom3  applies  to 
ions . 


The  plasma  temperature  can  be  estimated  by  determining  the  various 
species  in  the  railgun  plasma  and  then  using  the  method  of  PB  l 4)  to  evaluate 
the  degree  of  ionisation.  In  their  paper,  and  also  in  a  later  work  by  Powell 
151,  the  degree  of  ionisation  for  copper  is  calculated  as  a  function  of 
temperature  by  using  Saha's  equation  and  assuming  that  only  first  and  second 
ionisations  occur  in  the  plasma.  Using  this  facet  of  their  work,  the 
following  method  was  devised  to  estimate  the  plasma  temperature. 


If  the  atomic  masses  of  the  atoms  in  the  plasma  are  almost  the  same, 
then  momentum  transfers  1201  can  be  neglected  and  the  partition  of  energy 
between  the  different  elements  of  the  plasma  will  depend  largely  on  their 
respective  concentrations  assuming  each  of  the  respective  concentrations  of 
atoms  or  ion3  is  at  the  same  (average)  temperature.  Therefore,  spectral 
emissions  from  each  concentration  in  the  plasma  is  representative  of  not  only 
the  individual  species  but  the  plasma  as  a  whole. 


specif ically,  spectral  emissions  from  each  ion  concentration  win 
Indicate  the  degree  of  ionisation  reached  by  each  element  belonging  to  the 
plasma.  For  example,  we  know  from  PB’s  work  that  there  is  no  significant  2nd 
ionisation  below  16.3  x  io3  K  for  a  copper  plasma  at  10  MPa.  Therefore,  if  we 
only  detect  single  ionised  species,  then  the  peak  temperature  for  a  copper 
plasma  is  below  16.3  x  io3  K  at  10  MPa.  Conversely,  the  presence  of  2nd 
ionisation  means  that  the  peak  temperature  is  above  16.3  x  lo3  K. 

if  the  calculations  are  extended  to  include  elements  with  different 
ionisation  energies,  then  the  range  for  the  peak  plasma  temperature  is 
narrowed  by  identifying  the  presence  or  absence  of  ionised  species.  For  this 
reason,  a  number  of  elements  were  Introduced  into  the  railgun  plasma. 

These  elements  and  those  of  the  rails,  gun-body  and  the  atmosphere  provided  a 
range  of  ionisation  energies  from  5.21  ev(Ba+)  to  35.108  ev(o++).  For  each 
element  (ten  in  all),  the  degrees  of  first  and  second  ionisation  were 
calculated  at  a  fixed  pressure  using  Saha's  equation. 
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In  order  to  apply  Sana's  system  of  equations  to  a  rallgun  plasma  the 
following  assumptions  are  made: 

(1)  The  magnetic  field  does  not  greatly  influence  the  equations  and  the 
energy  levels  used  in  evaluating  the  partition  functions. 

(2)  A  rallgun  plasma  Is  in  "Local  Thermodynamic  Equilibrium"  ( LTE) , 
which  means  that  the  densities  in  specific  quantum  states  are  those 
pertaining  to  a  system  In  complete  thermodynamic  equilibrium  with 
the  same  total  mass-density,  temperature  and  chemical  composition  as 
the  actual  system.  LTE  can  be  expected  to  hold  for  dense  laboratory 
plasmas  where  colllsional  processes  involving  free  electrons  are 
predominant  1201. 

These  assumptions  are  discussed  In  conjunction  with  the  experimental  results 
in  section  4. 


where 


and 


Saha's  system  of  equations  151  is: 


_1±L 


x  .( 1 
j 


<v> 


5/2 


e 


-I  ,/k  T 
j  B 


K  (T,P) 
1  +  1 


(1) 


Xj  -  ratio  of  the  number  density  of  lon3  ionised  j  times  to  the 
total  number  density  of  heavy  particles  in  the  plasma, 
me  -  electron  mass, 

*B  -  Boltzmann's  constant, 

I.j  -  ionisation  energy  needed  to  ionise  the  atom  j  times 
p  «  the  plasma  pressure, 

o  -  ratio  of  number  density  of  electrons  to  the  number  of 
heavy  particles  =■  Ejx., 

T  -  the  absolute  temperature  of  the  plasma, 

Z-j  =  electronic  partition-function  for  the  ion  ionised  j-tlmes 
(i.e.  an  ion  in  the  j-th  ionisation  state). 


The  Xj's  are  normalised  so  that  their  sum  equals  unity,  i.e.  E  x.  -  1.  The 
electronic  partition  function  is  defined  by:  3  3 


3 


where 

and 

9ne 


-U  .,/lcT 
]i  B 


(2) 


-  energy  of  the  electronic  state  of  the  1-th  excited  level  of  an 
ion  in  the  j-th  ionisation  state 
is  the  degeneracy  factor  for  that  level. 


From  Equation  (1)  the  ionisation  of  the  j-th  ion  depends  on  the 
ionisation  of  the  (j+l)-th  Ion.  When  only  single  and  double  ionisation  is 
considered,  the  parameter  >  becomes: 


a 


+  2X. 


(3) 


So  that  Equations  (l)  become: 


10 


(4a) 


X  (x  +  2x2> 


(l-x1-x2>  (1+X1+2X2) 


-  K^T.P) 


and  X2 ( Xj +2X2 ) 

X, ( 1+x. +2x  ) 
1  12 

Equations  (4a)  and  (4b) 
given  in  PB: 


-  K2<T,P)  (4b) 

can  be  rewritten  in  the  following  forms  as 


2  K2  1  KjX2 


(5) 


and 


4K  ( 1-X  )  ,  ,, 


(2  +  3K,)x,  4K  ( l+K  ) ( 1+X  -2X  )  . 

1  2  (  { - ~ - “ - ^ -  +  li/2  -1,. 


l+K, 


X22(2+3Kl>2 


1  2 


(6) 


Equation  (6)  is  solved  iteratively  by  assuming  a  value  for  x2  on  the 
right  hand  side  and  then  calculating  a  new  x2  value.  The  process  is  repeated 
until  the  old  and  new  value  agree  to  the  desired  accuracy.  The  degree  of 
first  ionisation  Xj  follows  from  Equation  (5). 


PB  used  firstly  10,  and  then  15  terms  in  the  partition  function. 

They  found  no  discernible  differences  in  the  results.  However,  Grlem  120] 
states  that  all  levels  up  to  the  reduced  ionisation  level  have  to  be 
considered  because  for  higher  orbital  momenta  the  statistical  weights  are  high 
and  their  omission  will  lead  to  serious  errors  in  the  excited-state 
contribution  to  the  partition  function. 


To  avoid  any  potential  error,  three  sets  of  calculations  for  copper 
(the  computer  code  for  which  is  listed  in  the  Appendix)  were  carried  out  in 
which  the  fir3t  15  terms  were  used,  then  all  the  terms  up  to  the  ionisation 
level  and  finally,  all  the  documented  terms  in  Moore  (21).  The  iteration  was 
continued  until  successive  values  of  x2  agreed  to  within  1  x  10~5,  (as  in 
PB).  The  re3ult3  of  the  calculations  are  shown  in  Figure  6. 


The  graphs  for  the  relative  concentrations  of  first  and  second 
ionisation  obtained  by  truncating  the  sums  in  the  partition  functions  to  the 
first  15  energy  levels  matched  the  graphs  given  in  PB  and  hence  provided  a 
check:  on  our  calculations.  Incorporating  additional  excited  energy  levels 
modified  the  graphs  for  temperatures  above  13  x  lo3  K.  The  inclusion  of 
additional  energy  levels  up  to  the  ionisation  limit  of  each  species  lowered 
the  concentration  of  singly  ionised  atoms  and  raised  the  temperature  at  which 
the  pea*  in  first  ionisation  occurred.  However,  the  relative  concentration  of 
doubly  ionised  atoms  for  this  case  remained  unaffected.  The  inclusion  of  all 
the  energy  levels  documented  in  Moore  to  evaluate  the  partition  functions 
lowered  the  concentration  of  singly  ionised  atoms  for  temperatures  between 
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13  x  lo3  and  30  x  lo3  K.  For  temperatures  above  30  x  103  K,  the  concentration 
of  singly  ionised  atoms  was  above  that  given  in  PB  while  the  concentration  of 
doubly  ionised  copper  at  a  given  temperature  was  less. 


FIGURE  6  Level  of  ionisation  of  copper  versus  loge  (temperature) 
at  a  pressure  of  10  MPa. 


As  a  consequence  of  the  above  results,  further  calculations  in  this 
report  include  all  the  values  documented  by  Moore  (21)  even  though  Griem  1201 
states  that  only  levels  up  to  the  reduced  ionisation  energy  should  be 
considered.  Griem  was  concerned  that  in  summing  the  energy  levels  in  the 
partition  functions,  the  above  procedure  would  result  in  a  given  state  being 
summed  twice.  It  was  considered  that  as  the  levels  documented  by  Moore  are 
essentially  all  observed  values,  duplication  of  a  given  state  in  the 
summations  would  not  occur. 
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4. 


EXPERIMENTAL  RESULTS 


The  spectre  obtaihed  from  the  3  firings  are  shown  in  Figure  7.  Line 
reversals  and  the  copper  reference-spectrum  are  clearly  identified. 


A I 


Zn 


BaClj 

3 


Microphotometer  traces  made  from  the  films  (Figures  8a, 8b  and  8c) 
were  read  on  a  Farol-chart  reader.  Dispersion  (°A/cm)  of  the  wavelengths  on 
the  chart  records  was  determined  using  the  3273. 96°A  and  5218. 8°A  lines  in  the 
copper  reference  spectrum.  Razor  scratches  across  tne  reference  and  plasma 
spectra  provided  the  reference  line  from  which  readings  were  taken. 

Measurement  accuracy  was  ±o.5°A.  Line  positions  in  the  copper  reference- 
spectrum  were  determined  to  the  centre  line  at  the  half-height  level.  This 
method  was  also  used  for  measuring  the  position  of  the  line  reversals.  When 
it  was  not  possible  to  determine  the  half  height  (e.g.  for  small  line 
reversals),  the  peak  position  was  measured.  It  should  be  noted  that  Figures 
8a,  8b  and  8c  do  not  3how  the  detail  that  the  original  chart-recorder  records 
revealed. 


To  identify  an  ion,  the  measured  wavelengths  of  the  line  reversals 
were  compared  with  the  most  prominent  wavelengths  between  3000°A  and  6250°a 
for  the  excited  atom  or  ion  (Table  3).  These  lines  were  determined  from  line 
intensities  listed  in  References  117),  (23!  and  [241.  An  excited  atom  or  ion 
was  considered  to  be  present  if  the  three  most  prominent  lines  between  3000°A 
and  6250°A  were  identified.  This  is  considered  a  reasonable  approach  since 
the  transition  probabilities  117,221  for  the  Cull  and  Znll  lines  listed  in 
Table  3  have  a  similar  though  in  general  slightly  greater  magnitude  to  those 
for  Cul  and  Znl  respectively.  Assuming  that  line  absorption  occurs  for  the 
single-ionised  species,  it  follows  that  if,  for  example,  Znll  and  Cull  are  not 
detected  in  the  plasma,  then  their  concentration  are  below  the  concentrations 
of  the  detected  neutral  species.  A  similar  argument  is  expected  to  apply  to 
A1II  although  the  transition  probabilities  for  the  lines  considered  in  this 
case  are  not  available. 

The  wavelengths  measured  from  Firings  l,  2  and  3  are  tabulated  in 
Tables  4,  5  and  6  respectively  together  with  the  possible  identification  of 
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FIGURE  8c  Microphotometer  trace  of  Firing  No.  3:  Baci2  in  A1  foil 


the  lines.  Table  7  lists  the  identified  species  as  determined  by  the 
previously  stated  criteria. 

The  first  notable  feature  of  the  results  appearing  in  Table  7  is  the 
lacK  of  ionised  species,  which  was  unexpected.  Two  possible  explanations  for 
this  result  are: 


1.  The  plasma  was  a  partially  ionised  gas  composed  mainly  of 
excited  neutral  atoms. 

2.  The  absorption  spectra  obtained  in  these  firings  yielded  the 
lines  for  the  species  present  only  in  the  cooler  exterior 
regions  of  the  plasma  armature. 

If  the  first  explanation  is  valid,  then  a  low  temperature  estimate 
for  a  railgun  plasma  is  obtained,  as  shown  later  in  this  section.  The  second 
explanation  questions  the  use  of  absorption  spectra  to  determine  maximum 
temperature  estimates  in  an  optically  thlclc  railgun  plasma. 
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TABLE  3 


Prominent  lines  of  excited  elements  within  the  range  3000°A  to 
6250°A  using  References  [17,22,231. 


Excited  Atom* 


wavelengths  (°A) 

intensity  decreases  from  left  to  right 


cul 

Cull 

cum 

All 
All  I 

aiiii 

Znl 

Znll 

znm 

Bal 

Ball 

Balll 

01 

on 

om 


Nil 

Kill 

HI 

ci 

CII 

cm 

Cdl 

Cdll 

Cdlll 

CII 

cm 

cim 


3247.540,  3273.96,  S218.20,  3036.10,  3063.41 
4909.734,  3686.555,  5051.778,  4931.698 
3776.97,  4352.80,  4377.11,  3744.70 

3961.527,  3092.713,  3944.032 
3586.692,  3900.68,  3655.00, 

5696.60,  5722.73,  4529.19,  3601,63 

3345.020,  3302.588,  3302.94,  3282.33,  3345.57 
4911.664,  4924.03,  5894.33,  6021.18,  6111.53 
4970.8,  5579.0,  5563.8 

5535.551,  5777.65,  5519. OS,  3501.11,  5424.55 
4554.042,  4934.086,  5853.68,  4130.64 
3368.175,  3079.136,  3369.677,  3043.42 

6158.20,  6156.77,  6155.98,  5958.58 
4075.869,  4189.793,  4649.14,  3973.26 


3265.46, 

3759.87, 

3260.98, 

3961 . 59 

5725.50, 

4963 . 98, 

5829.54, 

5764 . 75  , 

3995.00, 

4630.54, 

5679.56, 

5005 . 14 

4379.11, 

4097.31, 

4103.4, 

4003.58,  5320.82 

4101 .73, 

3970.074 

,  3889.055,  3835.397 

5380.34, 

5052 .17, 

4932 .05, 

4771 . 75 

4267.26, 

4267.00, 

3920.677 

,  4076.00,  5145.16 

4647. 42, 

4650.25, 

5695.92, 

4651 . 47 

3610.51, 

3466.20, 

5085 . 82  , 

3403.65 

5378.3, 

5337.40, 

4415.63 

3035.72 

6140.25, 

6114.41, 

6194.757 

.  4526.19,  4438.49 

4794.54, 

5423.23, 

4810.06, 

4896.77,  5217.94 

3340.36, 

3191 .4, 

3139.16, 

3320.14 

*  Spectra  of  neutral,  single  and  double  ionised  atoms  are  indicated 
by  I,  II  and  III  respectively. 
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TABLE  4 


Measured  wavelengths  of  line  reversals  for  Firing  Ho.  l: 
Aluminium  foil  used.  (1  1  means  centre  line  at  half  height) 


wavelength  (°A) 

(±  O.S°A) 

Possible  Identification 

13247.7] 

3247.54  Cul 

13273.91 

3273.96  CUl;  3273.94  011 

13584.91 

3584 . 98  CII 

3603.3 

- 

3613.5 

3613.76  Cul 

3622 . 9 

- 

3634.0 

- 

13737.91 

- 

13947.71 

possibly  3944.0  All*;  3947.3,  3947.5,  3947.6  OI 

1 3954.91 

possibly  3961.5  All* 

4051.5 

4050.6  Cul 

14230 . 0 ] 

4230.4  Cull 

4677.6 

4678.156  Cdt 

4679.5 

- 

5081 . 8 

- 

5089 . 9 

5088.26,  5088.48,  5088.98  Cull 

5105 . 1 

5105.5  cm 

5106.7 

- 

5109.2 

5108.3  Cull 

[5106.8] 

- 

5129.9 

- 

5145.0 

5145.6  A1II,  5144.4,  5144.8,  5144.9  A1II;  5145.1  CII 

5218.8 

5218.2  CUI 

5784 . 0 

- 

15888.9] 

5889.97  CII 

5887.1 

- 

5888.4 

- 

5891 . 2 

5891.65  CII 

*  These  strong  lines  have  possibly  been  shifted  due  to  the  Stark  effect. 


17 


TABLE  5 


Measured  wavelengths  of  line  reversals  for  Firing  No.  2:  Zinc  foil 
used.  ((  1  means  centre  line  at  half  height) 


wavelength  (°A) 
(±0.5°A) 

Possible  Identification 

3056 . 8 

3057.1  All 

3063 . 7 

3063.4  Cul ;  3064.3  All 

3080 . 9 

- 

3091 . 7 

3093.9  cul ;  3092.713  All;  3092.710  A1II 

3100.3 

3099.9  CUI 

(3241.41 

3247.54  CUI 

(3274.31 

3273.96  CUI 

(3281.71 

3282.3  Znl;  3281.69  Cull;  3282.7  Cul 

(3303.2) 

3302.58  Znl;  3303.516  Cul 

(3346.41 

3345.5  Znl;  3345.9  Znl;  3345.02  Znl 

3360.0 

3361.09  CII 

3586.6 

3585.8  CII,  3586.69  A1II 

3721.7 

3720.77  CUI 

3724.2 

- 

388S.3 

3883.3  Znl;  3884.1  Cull 

3937.2 

3938.6  All I 

(3946.91 

3945.5  Cull;  3947.3,  3947.5,  3947.6,  OI;  3944.0*  All 

( 3964 . 8 1 

3965.4  Znl;  3961.5*  All 

3971  .  a 

- 

(4230.51 

4230.4  Cull;  4231.3  Cl 

4641.5 

4641.8  011 

4651 . 4 

4651.1  Cul;  4650.8  Oil;  4650.6  A1II 

4654.5 

4654.2,  4654.5  OI;  4655.05  A1II 

4675.3 

4674.76  Cull; 

( 4673 . 7 1 

- 

( 472S . 0 1 

- 

4737.7 

- 

(4813.11 

4812.9  Cull;  4812.8  Cl;  4810.5  Znl 

4846 . 1 

- 

4852.7 

- 

4868. 8 

- 

5080  4 

- 

These  strong  lines  have  possibly  been  shifted  due  to  the  Stark  effect. 
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Measured  wavelengths  of  line  reversals  for  Firing  No.  3:  Aluminium  wrapped 
around  Baci2.2H20  used,  (t  1  means  centre  line  at  half  height) 


wavelength  (°A) 

( ±0 . 5°A) 

Possible  Identification 

3020 . 4 

3057.3 

3057.154  All 

3065 . 0 

- 

3083.0 

3082.11  All;  3082.153  A1II 

3093 . 0 

3092.713  All;  3092.710  A1II 

13247.5) 

3247.54  CUI 

(3273.91 

3273.54  CUI 

(3302.1 1 

3302.588;  3302.94  Znl 

13345.14) 

3345.02;  3345.57  Znl 

(3502.2) 

3501.11  Bal 

3571.5 

- 

3581  .  7 

- 

3583 . 4 

- 

3736.5 

- 

(3948.0) 

3944.03  All;  3946.406*  A1II;  3947.3,  3947.5,  3947.6,01 

(3963.4) 

3961.5*  All 

(4230.5) 

4230.4  Cull 

4367.8 

4366.9  Oil;  4368.14  CII;  4368.3  OI 

(4463.6) 

- 

(4557.2) 

- 

4686.3 

4687.77  Cull 

4727.0 

4726.45  Bal;  4727.21  CII 

4814.3 

- 

4846.6 

- 

(4935.21 

4935.03  NI,  4934.086  Ball 

5082.0 

- 

5106 . 0 

5105.5  CUI 

5531.5 

- 

5881.9 

- 

5887.2 

- 

6101.5 

- 

6132.2 

- 

These  strong  lines  have  possibly  been  shifted  due  to  the  stark  effect. 


TABLE  7 


Atomic  species  detected  in  Firings  1-3 


Firing 

Foil  Material 

Atomic  Species  Identified 

1 

Aluminium 

Cul  and  possibly  All 

2 

Zinc 

Cul,  znl  and  possibly  All 

3 

Bacij  and  Aluminium 

Possibly  Cul,  All  and  Znl 

Another  feature  of  the  spectra  is  an  apparent  wavelength  shift  of 
the  very  strong  transition  lines,  3961. 5°A  and  3944. o°A  of  All.  For  example, 
in  Firing  No.  l,  aluminium  foil  was  used  and  it  would  be  reasonable  to  expect 
that  All  and/or  A1II  would  be  present.  There  were  two  strong  reversals  at 
3965. 9°A  and  3947. 7°A,  both  3llghtly  longer  than  the  3961. 5°A  and  3944. o°A 
lines.  These  two  lines  would  not  have  been  identified  but  for  the  fact  that 
in  Firing  No.  2  there  were  also  two  strong  reversals  near  the  same  region  at 
3946. 9°a  and  3964. 8°a.  In  addition,  when  aluminium  was  wrapped  around  the 
Baci2  in  Firing  No.  3,  there  were  2  strong  reversals  at  3946. o°A  and 
3963. 4°a.  The  lines  in  Firing  No.  2  if  attributed  to  All,  could  be  due  to 
contamination  of  the  railgun  bore.  The  mo3t  likely  source  of  aluminium 
contamination  would  be  aluminium  alloyed  or  welded  on  the  rail  surfaces  [211. 
Such  contamination  would  not  have  been  removed  by  the  3lmple  cleaning  method 
previously  described. 


The  reason  for  any  shift  in  these  two  wavelengths  may  be  due  to  the 
Stack  shift.  Such  shifting  of  the  centre  line  frequencies  is  due  to  charged 
particle  interaction  and  depends  on  the  atom  and  the  lines  involved  [26]. 
Because  of  the  repeated  occurrence  of  the  two  very  strong  line  reversals,  it 
is  likely  that  the  lines  are  due  to  All. 


Another  Interesting  result  was  that  Cdl  was  not  detected.  The 
absence  of  cadmium  can  be  explained  by  its  low  concentration,  as  cadmium  was 
only  0.6%  of  the  alloyed  rails.  Even  if  rail  material  constituted  half  of  the 
resultant  plasma,  there  would  only  be  0.3%  of  cadmium  within  the  plasma. 


Many  of  the  line  reversals  on  the  spectra  remain  unresolved  because 
they  could  not  be  Identified  with  any  of  the  species  expected  in  the  plasma. 
These  line  reversals  may  be  due  to  molecular  transitions.  Although  the 
unknown  line  reversals  were  checked  with  persistent  band  heads  for  a  large 
number  of  molecular  species  appearing  in  Pearce  and  Gaydon  1271,  not  one  was 
identified . 
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Having  Identified  some  of  the  species  which  were  present  within  the 
plasmas,  a  value  for  the  plasma  pressure  was  required  In  order  to  calculate 
the  temperature.  Since  no  pressure  measurements  were  made  during  each  firing, 
the  pressure  was  estimated  by  dividing  the  electromagnetic  force  on  the  plasma 
( L’ 1/2)  by  the  bore  cross-sectional  area,  as  given  In  McNab  13). 

Thus,  P  =»  ( I2L' ) / ( 2A)  (7) 


where  L'  »  Inductance  per  unit  length  of  the  rails  (H/m), 

I  -  current  through  the  rails  and  plasma  (A) 
and  A  =  cross-sectional  area  of  the  bore  (m2). 


It  is  shown  later  chat  the  predictions  for  the  degrees  of  Ionisation  using  the 
Saha  equations  are  not  significantly  affected  when  the  estimate  for  the 
pressure  of  a  railgun  plasma  obtained  from  Equation  (7)  is  Increased  or 
decreased  by  50  percent. 

The  equation  of  state  Is  given  in  PB  by  the  following  form:- 

V 

P«(l+a)p  -  ( 8 ) 

m 

o 

where  mQ  »  atomic  mass  of  the  Ion  or  neutral 

and  p  =  mass  density  of  the  arc. 

and  o  is  as  defined  in  Equation  (3). 


From  Equation  (8),  It  follows  that  the  maximum  temperature  occurs  when  the 
pressure  is  a  maximum,  provided  <l+o)  and  p  do  not  vary  substantially,  which 
Is  a  requirement  of  LTE  1201.  From  Equation  (7),  the  peak  pressure  occurs  at 
peak  current  and  hence  the  maximum  temperature  occurs  at  peak  current.  In 
Firing  No.  2,  the  peak  current  was  79.2  kA  as  3hown  In  Figure  9. 


TIME,  u* 

FIGURE  9  Current  record  for  Firing  No.  2 
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The  parameter  L'  in  Equation  < 7 )  was  found  following  Kerrislc's 
method  (28!  and  for  the  rails  used  In  this  series  Its  value  was  found  to  be 
0.41  ,H/m.  Therefore,  the  resultant  peaK  pressure  for  Firing  No.  2  was 
12.8  MPa . 


The  levels  of  the  degrees  of  first  and  second  ionisation  for  copper, 
zinc  and  aluminium  were  calculated  as  functions  of  temperature  using  a  plasma 
pressure  of  12.8  MPa.  The  results  are  shown  in  Figures  10,  ll  and  12 
respectively. 


Level  of  ionisation  for  copper  versus  loge<T) 
at  a  pressure  of  12.8  MPa 


FIGURE  10 


Because  the  level  of  sensitivity  necessary  to  detect  the  presence  of 
a  particular  species  could  not  be  determined,  we  consider  two  different  cases 
when  estimating  the  plasma  temperature.  For  the  first  case,  It  is  assumed 
that  lines  corresponding  to  concentrations  of  species  greater  than  0.5%  would 
appear  on  the  spectra.  For  the  second  case,  it  is  assumed  that  only  lines 
corresponding  to  species  concentrations  greater  than  10%  will  appear  on  the 
spectra . 


The  plasma  temperature  estimates  for  the  two  cases  of  0.5%  and  io% 
concentration  levels  are  presented  In  Table  8.  Single-ionised  aluminium  was 
not  found  in  the  plasma  armature,  which  Indicates  that  the  pealc  plasma 
temperature  was  less  than  6  x  lo3  and  ll  x  lo3  K  for  the  0.5%  and  10%  levels 
respectively.  Both  temperature  estimates  are  much  lower  than  the  theoretical 
estimates  given  In  References  (3,4,5  and  6).  These  low  estimates  for  the 
plasma  temperature  will  yield  low  values  for  the  conductivity  of  the  armature 
which  might  explain  why  the  theoretical  estimates  for  the  plasma  voltage  are 
lower  than  the  experimentally-determined  values.  The  discrepancy  between 
theoretical  and  experimental  values  for  the  plasma  armature  voltage  has  been 
attributed  to  large  voltage  drops  (-60  V)  at  each  rail-plasma  Interface  (31. 
If  the  plasma  temperature  is  low,  then  such  large  voltage  drops  may  not  be 
necessary  to  explain  the  measured  muzzle  voltage.  Recording  to  Reference 
(291,  the  cathode  and  anode  drops  are  only  of  the  order  of  10  V  for  copper 
electrodes . 


TABLE  8 

Temperature  for  the  onset  of  first  Ionisation  assuming  detection 
of  ionised  species  at  the  0.5%  and  10%  levels  of  ionisation 
for  a  pressure  of  12. a  MPa 


Plasma 

species 

Temperature 

<K> 

0.5%  level 

10%  level 

AIII 

6000 

1 1 ,000 

cull 

1000 

12,000 

znll 

7500 

0 

0 

0 

<M 

Because  the  rallgun  current  varied  during  each  firing  in  the  series 
and  the  value  of  the  rail  Inductance  gradient  (l')  was  not  accurately  Known, 
the  calculations  of  Saha's  system  of  equations  were  repeated  with  the 
pressures  deviating  ±  10%  and  ±  50%  (Figure  13)  from  the  original  pressure  of 
12.8  MPa.  The  results  are  shown  In  Table  9  for  only  the  10%  levels  of 
ionisation  for  cull  and  CuITI  because  the  0.5%  level  showed  no  discernible 
difference  In  temperature. 
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LEVEL  of  IONISATION 


NATURAL  LOO  of  TEMPERATURE 


figure  13  Level  of  ionisation  for  copper  versus  loge(T)  at  a 
pressure  of  6.4,  12.8  and  19.2  MPa 


TABLE  9 


Temperatures  for  Ionisation  levels  of  10*  in  copper 


Pressure 

(MPa) 

Temperature  (K) 

Cull 

cum 

6.4 

10,570 

23,700 

11.5 

11,170 

25,060 

12.8 

11,270 

25,310 

14.1 

11,360 

25,590 

19.2 

11,730 

26,450 
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As  Table  9  reveals,  a  +  50%  change  In  the  plasma  pressure  only 
produces  a  marginal  change  in  the  temperature  for  the  10%  level  of 
ionisation.  Figure  13  shows  that  at  a  pressure  of  6.4  MPa,  the  temperature 
ranges  over  which  first-  and  second-ionised  copper  are  expected  to  be  present 
are  almost  Identical  to  the  temperature  ranges  for  these  species  at  a  pressure 
of  12.8  MPa,  even  though  the  concentrations  vary  with  temperature.  The  same 
behaviour  is  found  to  apply  at  a  pressure  of  19.2  MPa.  Thus  changing  the 
pressure  by  a  factor  of  50%  does  not  have  a  significant  effect  on  ionised 
species  comprising  a  railgun  plasma.  Hence,  our  temperature  estimates  are  not 
very  sensitive  to  the  range  of  plasma  pressures  corresponding  to  variations  in 
railgun  current. 


4.1  Discussion  of  Assumptions 


The  major  assumption  used  in  determining  the  temperature  was  that 
the  plasma  was  in  LTE .  In  order  to  check  this  assumption  it  was  necessary  to 
estimate  the  plasma  particle  density.  Figure  14  shows  the  streak  photograph 
of  the  plasma  for  Firing  No .  3 .  In  this  figure  the  horizontal  dark  lines  do 
not  correspond  to  variations  in  plasma  intensity  but  are  due  to  position- 
markers  placed  along  the  gun-body.  Using  Figure  14,  the  plasma  armature 
length  was  found  to  be  35  ±  5  mm  at  peak  current.  As  unexploded  foil  was 
never  recovered  from  any  firings  at  MRL  [8,9,161,  this  suggests  that  nearly 
all  the  foil  became  plasma.  Even  if  only  10%  of  the  foil  became  plasma,  a 
number  density  greater  than  1024  m-3  is  obtained  for  all  the  firings.  This 
density  is  typical  for  dense  plasmas  [311  and  indicates  that  colllslonal 
processes  involving  free  electrons  may  dominate,  which  Is,  as  mentioned 
previously,  a  criterion  for  the  validity  of  LTE. 


FIGURE  14 


400  *3  streak  photograph  of  Firing  No  3.  Graticule 
divisions  are  8  »s  apart. 
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A  second  assumption  used  in  obtaining  the  temperature  was  that  the 
effect  on  the  spectrum  of  the  strong  magnetic  field  In  the  plasma  armature  was 
negligible.  To  show  this,  the  spectral  trace  for  Firing  3  was  examined  for 
evidence  of  the  Zeeman  effect.  According  to  Griem  (321,  the  Zeeman  effect 
should  produce  a  wavelength-shift  for  Gaussian  profiles  of:- 

-8  2 

AX  =  10  X  B 


where  x  Is  the  wavelength  (°A)  and  B  is  the  magnetic  field  (T). 

The  maximum  magnetic  field  obtained  by  assuming  a  thin  current  sheet 
distribution  in  the  rails  (30]  is  found  to  be  about  6.<  T.  The  Zeeman  shift 
for  the  All  line  at  3961.4  °A  is  therefore: 


AX  =  l°A 


Since  the  expected  Zeeman  shift  is  very  small,  magnetic  field  effects  are 
unlikely  to  interfere  with  the  validity  of  Saha's  equations  for  a  railgun 
plasma . 


It  should  also  be  mentioned  that  time  variations  and  spatial 
inhomogeneities  in  the  plasma  may  further  restrict  the  validity  of  lte.  in 
order  to  investigate  these,  the  streak  film  was  colour-enhanced  by  using  an 
image-analysis  system  developed  at  MRL  (32).  The  intensity  contours  resulting 
from  this  analysis  are  shown  in  Figure  15. 


figure  15  Colour-enhanced  image  showing  intensity  contours  of  the  plasma  in 
Firing  Ho.  2.  The  intensity  contours  surrounding  the  plasma 
region  before  t-o  are  believed  to  be  due  to  flaring  in  the  image 
lntensifler  tube. 
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As  shown  in  Figure  15  intensity  gradients  appear  at  the  rear  and 
leading  edge  of  the  expanding  plasma.  These  intensity  gradients  show  the 
cool  layers  of  plasma  which  are  predicted  by  theory  1 4).  The  intense  central 
region  (green)  reveals  a  hotter  region  of  the  plasma,  which  is  stable  in 
Intensity  for  approximately  30  *s  around  peak  current  or  peak  temperature.  It 
is  expected  that  any  significant  variation  in  the  internal  temperature  and 
density  of  the  plasma  would  be  indicated  by  intensity  gradients  appearing  in 
the  central  region,  which  did  not  occur.  Thus  photographic  evidence 
indicates  that  a  hot  stable  region  exists  within  the  plasma,  thereby 
supporting  additional  criteria  for  the  validity  of  LTE. 


5.  ANALYSIS  OF  PREVIOUS  TEMPERATURE  ESTIMATES 


The  temperature  estimate  obtained  in  the  previous  section  by 
applying  Saha's  eguation  to  the  various  species  identified  from  absorption 
spectra  is  much  lower  than  the  peak  temperature  estimates  found  in  previous 
theoretical  and  experimental  work  except  for  the  work  of  Gathers  and 
Hord  (101.  Unfortunately,  Gathers  and  Hord  did  not  report  the  circuit 
parameters  or  railgun  si2e,  so  it  is  difficult  to  compare  their  result  with 
the  results  in  this  report. 


5.1  Richardson's  Emission  Spectroscopy 


In  order  to  establish  whether  absorption  spectra  do  provide  a  valid 
estimate  for  the  peak  temperature  of  a  railgun  plasma,  our  method  was  applied 
to  the  ions  and  excited  neutral  atoms  identified  from  the  emission  spectra 
obtained  by  Richardson  (8).  Richardson's  railgun  had  a  smaller  bore  (48  mm2 
compared  with  loo  mm2)  but  the  peak  plasma  current  was  almost  identical  (80 
kA) .  Since  our  method  for  determining  plasma  temperature  differs  from 
Richardson's  approach,  it  was  considered  of  interest  to  see  if  the  two 
different  approaches  yielded  similar  temperature  estimates. 


In  carrying  out  this  comparison  the  same  rule  for  identifying  the 
neutral  and  ionic  species  in  the  absorption  spectra  was  applied  to 
Richardson's  results  (8).  The  three  most  prominent  line-emissions  of  a  given 
ionic  species  within  the  observed  spectral  range  were  sought  (Table  lo).  It 
was  found  that  Cul,  All  and  Oil  and  possibly  Cull,  Alii  and  Nil  were  present 
in  the  muzzle  flash.  However,  OI  and  Nl  were  not  found. 


The  approach  used  by  Richardson  to  obtain  a  peak  temperature 
estimate  was  to  select  the  ionic  species  with  the  highest  ionisation  potential 
in  the  muzzle  flash  and  then  to  estimate  the  plasma  temperature  by  evaluating 
the  temperature  at  one-fifth  of  the  ionisation  energy  for  that  species  as 
Parkes  and  Strachan  (2)  had  done.  Richardson  chose  Oil  and  hence  obtained  a 
temperature  greater  than  30  x  103  K.  We  believe  that  this  approach  is 
inappropriate  because  it  neglects  the  range  of  temperatures  over  which  ionic 
species  can  exist  in  a  plasma  (Fig  16).  In  addition,  Fig.  16  reveals  that  for 
a  temperature  of  30  x  103  K,  at  least  80  percent  of  the  oxygen  in  the  plasma 
would  be  doubly  ionised.  Lines  corresponding  to  OIII  were  not  found. 
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TABLE  10 


The  Most  Prominent 

spectral  Lines  for  Different  Species  117) 
within  the  spectral  range  <3360  to  4750°A)  observed  by  Richardson  181 


Excited  Atom 

intensity 

wavelength  (°A) 
decreases  from  left  to  right 

All 

3961.52, 

3944.0,  3443.64, 

All  I 

3900 . 68 , 

3587.07,  3655.0 

Aim 

4529.19, 

3600.6,  4512.5 

cm 

3530.38, 

4062.64,  3599.13, 

Cull 

4043.48, 

4073.7,4161.14 

cum 

3776 . 97, 

4352.80,  4377.11 

01 

3947.29, 

3947.48,  3947.59 

on 

3973.26, 

4075.86,  4189.79,  3911.96 

OIII 

3759.89, 

3961.59,  3754.65 

NATURAL  LOG  of  TEMPERATURE 

figure  ic  ionisation  of  oxygen  versus  loge<T)  at  a  pressure  of  12.8  MPa 


The  graphs  for  the  first  and  second  degrees  of  Ionisation  of  oxygen 
using  all  the  energy  levels  documented  In  Moore  121)  appear  In  Figure  16. 
According  to  this  figure,  the  degree  of  first  ionisation  peaks  at  about  90% 
when  the  temperature  Is  about  36  x  io3  K.  At  this  temperature  Cull,  Cum, 
Alii  and  A1III  would  be  expected  to  be  present  In  the  plasma  as  Indicated  by 
Figures  10  and  12.  Since  the  relative  concentrations  of  01  and  OIII  are  small 
at  36  x  lo3  K,  it  is  feasible  that  spectral  lines  corresponding  to  oi  and  om 
cannot  be  detected.  On  the  other  hand,  according  to  Figure  10,  the 
concentrations  of  cull  and  CuIXI  are  35%  and  65%  respectively  at  36  x  103  K. 
However,  CuIII  was  not  found  even  though  a  copper  foil  was  used  to  initiate 
the  plasma  and  there  was  doubt  whether  Cull  was  detected.  These 
considerations  suggest  that  the  plasma  could  not  have  reached  a  temperature  of 
36  x  lo3  K.  Hence,  determining  the  temperature  from  Oil  alone  cannot  be 
considered  reliable.  These  conclusions  are  further  substantiated  by  the 
absence  of  A1III  and  possible  absence  of  A1II  lines  on  the  spectra*  since  A1II 
and  A1III  have  similar  concentrations  to  Cull  and  CuIII  at  36  x  lo3  k. 

Disregarding  Oil  yields  different  temperature  estimates  depending  on 
whether  cull  and  Alii  Ions  are  present  In  the  plasma,  if  these  ions  are  not 
present,  then  the  species  identified  in  the  emission  spectra  are  consistent 
with  those  in  our  absorption  spectra  and  the  maximum  temperature  estimate  of 

11  x  lo3  K  obtained  in  the  previous  section  are  valid  for  Richardson's 
experiments.  However,  if  both  Cull  and  Alii  ions  were  in  the  plasma,  then  the 
temperature  estimates  are  modified.  Assuming  that  the  concentrations  of  Cull 
and  CuIII  must  be  greater  than  0.5%  before  being  detected,  the  plasma 
temperature  estimate  ranges  from  7  x  103  to  20  x  lo3  K  according  to  Figure 

10.  Assuming  also  that  the  concentrations  of  A1II  and  A1III  ions  must  be 
greater  than  o.5%  before  being  detected,  the  temperature  according  to  Figure 

12  ranges  from  6  x  lo3  to  19  x  103  K.  Since  there  i3  some  evidence  that  both 
A1II  and  Cull  may  be  present,  this  restricts  the  temperature  range  to  between 
7  x  lo3  and  19  x  io3  K.  On  the  other  hand,  if  the  concentrations  of  the  four 
species  must  be  greater  than  10%  before  being  detected,  then  a  different 
plasma  temperature  estimate  is  obtained.  According  to  Figure  10,  the 
temperature  for  Cul  and  Cull  to  be  present  at  concentrations  greater  than  10% 
In  the  muzzle  flash  of  the  plasma  ranges  from  11  x  103  to  26  x  io3  k.  From 
Figure  12,  the  temperature  for  All  and  A1II  to  be  present  at  greater  than  this 
concentration  ranges  from  10  x  io3  to  25  x  io3  K.  Hence  the  common 
temperature-range  is  from  11  x  103  to  25  x  103  K,  which  is  higher  than  the 
range  obtained  by  assuming  that  the  concentrations  of  the  species  had  to  be 
greater  than  0.5%.  Since  emission  spectroscopy  is  very  sensitive,  the 
concentrations  of  A1II  and  Cull  are  almost  certainly  small.  Then  the 
temperature  estimate  is  closer  to  11  x  io3  K  rather  than  25  x  io3  k.  The 
temperature  estimates  presented  in  this  paragraph  are  below  all  previous 
theoretical  predictions  and  Richardson’s  minimum  peak  temperature  of 

30  X  103  K. 

Our  temperature  estimates  using  Richardson's  emission  spectra  were 
determined  at  the  pressure  of  12.8  MPa  (peak  current  of  80  kA)  whereas  his 


*  Aluminium  species  were  expected  in  the  plasma  because  Richardson  had 
Inserted  aluminium  plugs  into  the  copper  rails. 


observations  were  at  the  muzzle  when  the  pressure  (current)  was  lower. 

However,  It  is  shown  in  Section  3  that  large  pressure  changes  produce  marginal 
changes  at  the  ionisation  levels  considered.  Furthermore,  Fig.  13  shows  that 
the  temperature  decreases  as  the  pressure  decreases. 


S.2  Marshall's  Temperature  Estimate 

As  described  briefly  in  section  1,  Marshall  (9)  has  used  a  different 
approach  to  obtain  an  estimate  of  the  plasma  temperature.  However,  his  method 
does  not  consider  the  following  points: 

(1)  Plasma  flow  through  a  nozzle  is  affected  by  viscosity  which  is 
temperature  dependent,  varying  between  T1  2  for  a  lightly-ionised 
plasma  and  T5'2  for  a  highly-ionised  plasma  (331. 

(2)  The  "degrees  of  freedom"  associated  with  the  change  in  the  number  of 
particles  (dissociation,  association,  ionisation,  chemical  reactions) 
change  dramatically  for  an  expanding  and  cooling 

plasma  (341.  This  affects  the  ratio  of  specific  heat3,  which  appears 
in  Marshall '3  equation. 

Consequently,  we  believe  that  more  investigation  is  required  before 
this  method  can  be  considered  as  a  reliable  means  of  estimating  the  plasma 
temperature . 


«.  COMDOCTXVITY  CAUTOATIOIS 

From  the  absorption  spectra,  a  pealc  plasma  temperature  less  than 
11  x  103  K  (at  the  10%  sensitivity  level)  was  obtained,  whereas  from  the 
emission  spectra  a  pea*  temperature  as  high  as  25  x  103  K  was  obtained. 
Another  method  of  estimating  the  pealc  temperature  range  is  to  calculate  the 
electrical  conductivity  of  the  plasma  for  varying  temperatures  and  then  to 
compare  it  with  conductivity-values  obtained  by  using  muzzle  voltage  records 
and  Image-analysis  results. 


The  electrical  conductivity  of  a  raligun  plasma  was  calculated  by  PB 
using  a  modified  form  of  an  expression  of  Spltzer  1351.  The  expression  given 
In  PB. 


However,  this  approach  was  not  used  here  for  a  number  of  reasons. 
Firstly,  In  deriving  this  equation  It  Is  assumed  that  the  gas  Is  fully 
ionised.  In  addition,  Powell  and  Batteh  have  assumed  that  the  plasma  is 
nearly  completely  double-ionised  and  thus  have  set  Z  equal  to  a  value  of  2. 
This  may  not  necessarily  be  the  case  for  a  raligun  plasma.  Also,  according  to 
Cohen  et  al  (361,  the  theory  used  to  obtain  the  expression  used  for  plasma 
conductivity  by  PB  breads  down  completely  for  plasmas  with  high  densities  and 
low  Kinetic  temperatures.  In  view  of  the  doubtful  applicability  of  spitzer's 
expression  to  raligun  plasmas,  a  computer  code  developed  by  Kovitya  ( 3 v , 3 8  I 
wa3  used  Instead.  This  code  calculates  various  properties  of  partially  and 
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fully-lonlsed  plasmas  including  the  electrical  conductivity  and  has  been  used 
by  Kovltya  and  his  co-workers  In  their  studies  of  ablation-dominated  arc3  and 
sun-spot  activity.  They  report  good  agreement  with  experiment  In  References 
1  39-41 1 . 


The  electrical  conductivity  for  a  copper  plasma  (Fig.  n)  was 
calculated  using  the  Kovltya  code  for  temperatures  between  8  x  lo3  and 
26  x  io3  K  and  for  pressures  between  4  and  14  MPa.  Only  a  copper  plasma  was 
considered  because  the  degrees  of  Ionisation  for  copper,  aluminium  and  zinc 
with  respect  to  temperature  are  very  similar  (Figures  10,  11  and  12). 
Furthermore,  in  the  RAPID  Plasma  Intensity  Profiles  series  of  rallgun  firings 
114],  It  was  found  that  the  muzzle  voltage  records  were  very  similar 
regardless  of  whether  a  copper,  zinc  or  aluminium  foil  was  used  to  generate 
the  plasma.  The  pressure  range  Included  the  plasma  pressures  previously 
considered  in  Section  5.  The  graph  shows  that  for  temperatures  lower  than 
IS  x  103  K,  the  plasma  conductivities  are  comparable  whereas  for  temperatures 
above  IS  x  103  K,  the  plasma  becomes  more  conductive  as  pressure  increases. 


TEMPERATURE,  K 


figure  17  Electrical  conductivity  of  a  copper  plasma  vs  temperature 
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To  calculate  the  conductivity  of  the  plasma  from  experimental 
results,  the  voltage  across  the  plasma  at  peaX  current  was  obtained  from  the 
muzzle  voltage  record  shown  In  Figure  18. 


TIME,  fj  s 

figure  1#  Muzzle  voltage  versus  time  (Firing  No.  4) 


The  average  muzzle  voltage  at  peaX  current  for  a  number  of  firings 
was  230  V.  Allowing  a  total  electrode  voltage  drop  of  20  v  yields  a  plasma 
voltage  of  210  v,  which  when  divided  by  the  peaX  current  of  79  XA  produces  a 
plasma  resistance  of  2.7  mQ.  Since  the  plasma  length  was  found  to  be  75  ±  5 
mm  (Section  4.1)  and  the  area  of  the  bore  Is  100  mm,  the  value  for  the 
electrical  conductivity  of  the  plasma  is  (5.0  ±  0.3)  x  103  s/m.  From  Figure 
17,  the  corresponding  peaX  plasma  temperature  for  a  pressure  of  12.8  MPa 
occurs  Is  (10.7  t  0.3)  x  103  K,  which  falls  In  the  temperature  ranges 
estimated  from  our  spectroscopic  observations.  Even  If  each  electrode  voltage 
drop  Is  assumed  to  be  «ov,  a  peaX  temperature  of  only  14.7  x  103  k  (at  12.8 
MPa)  Is  obtained.  Thus  the  temperatures  found  by  this  method  are  closer  to 
the  temperature  estimates  obtained  from  the  absorption  spectra  than  those 
obtained  from  Richardson's  observations. 


It  should  be  noted  that  the  method  employed  In  this  section  was  used 
by  KowalenXo  (14)  to  obtain  a  plasma  temperature  of  20  x  103  K  for  the  RM- 
experlment  In  which  the  current  was  300  XA.  This  suggests  that  the 
temperature  of  a  rallgun  plasma  Is  not  a  strong  function  of  the  circuit 
current,  it  should  also  be  noted  that  if  the  actual  length  of  the  plasma  was 
only  one  half  of  the  observed  length,  then  the  plasma  temperature  would  be 
less  than  14  x  103  K,  which  Is  a  similar  result  to  KowalenXo  (14). 
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Furthermore,  the  lower  temperature  and  subsequent  higher  Lisistlvlty  of  the 
plasma  Indicates  that  the  assumed  plasma  rail  Interface  voltages  13,4)  are  not 
necessarily  so  high  and  much  lower  values  are  possible. 


7.  COMCLUSIOW 


In  this  report  Saha's  system  of  equations  has  been  used  to  find  the 
relative  concentrations  of  singly  and  doubly  ionised  atoms  as  a  function  of 
temperature  for  a  number  of  elements  expected  in  railgun  plasma  armatures. 

When  additional  energy  levels  were  used  to  evaluate  the  partition  functions  in 
Saha's  system  of  equations,  the  concentrations  of  Cull  and  CuIII  differed 
slightly  from  those  given  in  PB. 


The  first  method  of  determining  the  plasma  temperature  presented  in 
this  report  depends  on  the  identification  of  species  comprising  a  railgun 
plasma  appearing  in  absorption  and/or  emission  spectra.  This  method  is 
limited  by  the  detection  sensitivity  and  by  the  documentation  of  excited 
energy-levels  for  various  species. 


Only  excited  states  of  copper,  zinc  and  aluminium  were  identified  on 
the  absorption  spectra  of  the  free-flowing  railgun  plasma.  If  the  sensitivity 
level  of  the  absorption  spectra  was  such  that  only  concentrations  of  each 
species  greater  than  o.5»  would  be  detected  by  the  spectroscope,  then  the 
temperature  of  the  plasma  armature  was  found  to  be  les3  than  6  x  103  K.  if 
the  sensitivity  level  required  concentrations  greater  than  lo%  (the  more 
likely  case)  then  the  temperature  was  found  to  be  less  than  ll  x  103  K. 

Two  possible  explanations  were  given  for  the  absence  of  ionised 
species  in  the  plasma  and  hence  the  low  temperature  estimates.  The  first  was 
that  the  plasma  was  only  partially-ionised  and  thus  was  composed  mainly  of 
excited  neutral  species.  The  second  explanation  was  that  the  absorption 
spectra  contained  the  line  reversals  which  reveal  only  conditions  in  the 
cooler  regions  of  the  plasma.  In  an  attempt  to  determine  which  of  the  two 
explanations  was  valid,  the  emission  3pectra  of  the  muzzle  flash  obtained  by 
Richardson  18]  were  examined.  It  was  shown  that  the  method  used  by  Richardson 
of  estimating  plasma  temperatures  based  solely  on  the  presence  of  Oil  in  the 
emission  spectra  may  not  be  reliable.  If  there  were  few  singly-ionised 
species  present  in  Richardson’ 3  plasma,  then  the  temperature  estimates 
obtained  from  the  absorption  spectra  reported  here  are  consistent  with  the 
estimates  obtained  from  Richardson's  emission  spectra.  However,  some  lines  on 
the  Richardson's  emission  spectra  suggested  the  possibility  of  some  singly- 
ionised  species  being  present.  Allowing  for  these,  the  temperature  was  found 
to  range  from  7  x  lo3  to  19  x  lo3  K  (most  probable)  or  from  ll  x  lo3  to  25  x 
103  K  depending  on  the  assumed  detection  sensitivity. 


In  order  to  narrow  the  temperature  estimates  obtained  from  the 
spectral  analysis,  the  electrical  conductivity  was  calculated  with  respect  to 
temperature  and  then  compared  with  an  experimentally-determined  value.  The 
theoretical  values  for  the  electrical  conductivity  were  obtained  from 
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Kovltya's  computer  code  137,381.  This  code  was  used  because  the  Spltzer 
expression  for  the  electrical  conductivity  was  of  doubtful  applicability  to 
railgun  plasmas  and  because  the  code's  results  have  been  shown  to  give  good 
agreement  with  experiment  139-411.  comparison  of  the  theoretical  and 
experimentally-determined  values  for  the  electrical  conductivity  produced  a 
peak  plasma  temperature  of  (10.7  ±  0.3)  x  103  K,  which  was  in  agreement  with 
the  estimate  obtained  from  the  absorption  spectra  and  within  the  ranges 
obtained  from  the  analysis  of  the  emission  spectra.  This  suggests  that 
absorption  spectra  can  be  used  to  obtain  an  estimate  for  the  plasma 
temperature . 


In  conclusion,  our  examination  of  the  spectra  obtained  from  railgun 
firings  and  comparison  of  the  electrical  conductivity  with  results  obtained 
from  muzzle  voltage  records  and  image  analysis  indicate  that  the  temperature 
of  a  railgun  plasma  is  less  than  11  x  103  K  which  is  much  lower  than  other 
estimates  except  for  those  given  in  References  (10  and  14).  Substantial 
variation  in  the  assumed  pressure  did  not  affect  the  temperature  estimates 
greatly.  The  low  temperature  estimates  obtained  in  this  report  indicate  that 
the  electrical  conductivity  and  degree  of  ionisation  of  a  railgun  plasma  are 
much  lower  than  the  values  given  in  PB.  The  subsequent  higher  plasma 
resistivity  also  mean3  that  the  assumed  plasma/rail  interface  voltages  13,4) 
may  not  be  as  high  as  previously  thought. 
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APPENDIX  A 


PROGRAM  ION 

C  THIS  PROGRAM  CALCULATES  THE  LEVEL  OF  IONISATION  IN  A  PLASMA. 

C  THE  PLASMA  IS  CONSIDERED  TO  BE  IN  LOCAL  THERMODYNAMIC 

C  EQUILIBRIUM  (LTE).  THE  1ST  AND  2ND  DEGREES  OF  IONISATION 

C  ARE  CALCULATED  FROM  A  METHOD  DETAILED  IN  A  PAPER  BY  J.  D.  POWELL 

C  AND  J.  H.  BATTEH  CALLED  PLASMA  DYNAMICS  OF  THE  ARC  DRIVEN 

C  RAILGUN,  SEPT  1980. 

C 

C  THIS  PROGRAM  READS  A  GIVEN  DATA  FILE  AND  OUTPUTS  A  FILE 

C  WHICH  CAN  BE  READ  BY  THE  PLOTTING  ROUTINE  WRITTEN  BY  STEVEN 

C  KENNETT.  THE  1ST  AND  2ND  DEGREES  OF  IONISATION  ARE  PLOTTED 

C  (FOR  A  CONSTANT  PRESSURE)  AGAINST  THE  NATURAL  LOG  OF  THE 

C  TEMPERATURE . 

C  PROGRAM  RUN  ON  A  VAX  11/780  USING  VAX  FORTRAN  77. 

C 

C  THIS  PROGRAM  WAS  WRITTEN  BY  GREGORY  CLARK. 

C  MATERIALS  RESEARCH  LABORATORIES 

C  JULY,  1985. 

C  PROGRAM  SWITCHES  ....USE  /G  FLOATING  SWITCH 

£********* **************************************************** ****** 
REAL*8  Z(  3,30  0) ,T(200) , PRESS , KSAHAt 2 ) ,PART<3> , IP  <  2  > 

REAL* 8  U(3,300) ,DG( 3,300 ),FN( 6) ,Xl( 200) ,X2( 200) , CHECK 
REAL*8  CPI,C,A1,A2,A3,A4 ,QN, TEMP , B1 , B2 , B3 
REAL  ATOM(6),ANS 
INTEGER  NLEVEL ( 3 ) 

Q* ************************  ******  ************************************ 

C  THE  CONSTANTS  USED  IN  THE  PROGRAM  ARE  LISTED  BELOW  IN  THE 

C  FOLLOWING  ORDER:  CME-ELECTRON  MASS  ( KG ) , CH-PLANCKS 

C  CONSTANT  (JOULES  SEC ) , CKB-BOLTZMANNS  CONSTANT  (JOULES/KELVIN), 

C  CPI-PI 

C 

C  REFERENCE  USED .  FAYE  C  LABY . 14TH  EDITION 

Q  **  *  *  *  **************  *****  **  it****  **  ****  ****  ***  *  *  ****  **  **  *  *  ****  ***  * 

CME-9. 109534D-31 
CH-6 . 626176D— 34 
CKB-1 . 30O662D-23 
CPI-3.14159 

******************************* ***************** *********** ****** 

C  THE  FOLLOWING  CALCULATION  IS  A  CONSTANT  IN  THE  SAHA  EQUATION. 

C  IT  IS  CALCULATED  IN  THESE  STEPS  TO  REMAIN  WITHIN  THE  NUMBERS 

C  LIMIT  OF  THE  VAX  COMPUTER . 

Q** ***************************** *************** ********** ***** ****** 

C-2*CPI*(CME/CH)*(CKB/CH) 

C-C**l . 5 
C-2*C*CKB 

Q******************** ****** **** **** ***** **** *****"*** **************** 

1  TYPE  *,'  WHICH  DATA  FILE  IS  REQUIRED?’ 

ACCEPT  2 , FN 

2  FORMAT! 6A4) 

C  INPUT  PRESSURE  IN  PASCALS  (DOUBLE  PRECISION  FORMAT) 

TYPE  *,’  WHAT  PRESSURE?’ 


ACCEPT  *, PRESS 

OPEN( UNIT-1 , NAME-FN , TYPE- ' OLD ' , READONLY) 

OPEN(UNIT-2,TYPE-'NEW ,NAME-' ION. DAT’ ) 

qXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C  READ  NAME  OF  THE  ELEMENT 

READ (1,2) ATOM 

C  READ  IONISATION  POTENTIALS  (1ST.  t  2ND  IN  eV) 

READ(l,*)lP(l),IP(2) 

C  READ  THE  ATOMIC  DATA  FOR  THE  THREE  LEVELS  OF  IONISATION 

DO  30  M-1,3 

C  NOW  READ  THE  VALUE  OF  THE  QUANTUM  NUMBER  J  AND  THE  ASSOCIATED 

C  ENERGY  LEVEL  AND  THEN  CALCULATE  THE  DEGENERACY  FACTOR. 

DO  20  1-1,300 
READ(1,*)QN,U(M,I) 

C  DETECT  END  OF  DATA  FOR  A  GIVEN  IONISATION  LEVEL 

IFfQN.EQ. 1000. ) GOTO  21 

C  CONVERT  ENERGY  VALUES  (WHICH  ARE  DIRECT  FROM  C.E.  MOORE) 

C  INTO  eV  AND  THEN  INTO  JOULES . . 

U(M,l)-U<M,l)*l -2395D-«*1 . 60219D-19 
DG  ( M ,  I )  —  (  2  .  *QN )  + 1 . 

20  CONTINUE 

21  NLEVEL(M)-I-1 

30  CONTINUE 

C  CONVERT  eV  TO  JOULES  (IONISATION  POTENTIAL) 

IP(l)-IP(l)*l. 60219D-19 
IP(2)=IP(2)*1.60219D-19 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C  NOW  WE  CALCULATE  THE  1ST  AND  2ND  DEGREES  OF  IONISATION  FOR 

C  A  TEMPERATURE  RANGE  8,000K  TO  (NUM*500 1+7500K  IN  STEPS  OF  500K 

0XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

NUM— 1 70 

DO  1000  IT-1, NUM, 1 

TEMP— ( DBLE( IT) *500 . 1+4500. 

C****** ****************************************************** ****** 

C  NOW  WE  CALCULATE  THE  PARTI ION  FUNCTION  FOR  A  GIVEN 

C  TEMPERATURE . 

C 

C** ********************************************** ****************** 
Z(l,l)-0.0 
Z( 2, 1 )-0 . 0 
Z( 3 , 1 )— 0 . 0 
DO  100  JJ-1,3 
DO  100  1-1 , NLEVELf  JJ ) 

Z(JJ,I+1)-Z(JJ,I)+(DG( JJ, l)*DEXP(-U( JJ,I)/(CKB*TEMP) ) ) 

100  CONTINUE 

PART ( 1 ) -Z ( 1 , NLEVEL ( 1 ) + 1 ! 

PART ( 2 ) -Z ( 1 , NLEVEL ( 2 ) + 1 ) 

PART ( 3 ) -Z( 3 , NLEVEL ( 3 ) +1 > 

£************************** **v^  ************************************* 

C  CALCULATE  THE  VALUE  OF  SAHA'S  EQUATION. 

£XXXXXXXXXXXXXXXXXXX**XXXXXXXXXX*XXXXXX***XX*XXXX*X**XXXXXXXXXXXXXX 

DO  150  N-1,2 

KSAHA( N )-( C/ PRESS ) * ( PART ( N+l ) /PART(N) ) * ( TEMP** 2 . 5 ! 

ksaha(n)-ksaha(n)*dexp(-ip(n)/(ckb*temp) ) 

150  CONTINUE 

c* ***************************************************************** 
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C  CALCULATE  BY  ITERATION  THE  2ND  DEGREE  OF  IONISATION. 

C  LET  THE  INITIAL  VALUE  OF  IONISATION  (X2(lT)  )  EQUAL  0.01. 

C  4  00  ITERATIONS  HAVE  BEEN  ALLOWED  FOR.  THE  ITERATION  WILL 

C  END  IF  SUCCESSIVE  VALUES  OF  THE  2ND  DEGREE  DIFFER  BY  LESS 

C  THAN  1  PART  IN  100,000. 

C 

C****************************************************************** 

X2  ( IT ) =0 .01 
DO  300  L-1,400 
CHECK«X2<IT> 

Al-4  . *KSAHA<  2 ) *< 1 . -X2 (IT)) 

A1=A1/ ( KSAHA( 1 ) *X2 ( IT ) ) 

Al-(l.+Al)**0.5 
Al-Al-1 . 0 

£******************************************★**■************** 

C  THE  FOLLOWING  "IF"  STATEMENT  WAS  INSERTED  TO 

C  COVER  THE  INSTANT  WHEN  Al-1 . 0  GOES  SMALLER  THAN 

C  THE  15  DIGIT  LIMIT  WITH  THE  DOUBLE  PRECISION 

C  /G  SWITH. 

£***** ***** **************************************** ****** *** 
IF(Al.EQ.O.O)  TYPE  *,'A1  EQUALLED  0.0' 

IF(Al.EQ.O.O)  Al-l.D-16 
A1-(KSAHA(  1  )/KSAHA(2)  )*A1 
A2-2.+(3.*KSAHA(l)) 

A2-A2*X2(IT)/<1.+KSARA(1>> 

A3-4 . *KSAHA< 1 ) *( 1 . +KSAHA( 1 > ) 

A3-A3*(l.+X2(IT)-(2.*X2(IT)**2. )) 
A3-A3/(X2(IT)*(2.+(3.*KSAHA(1> ) ) >**2. 

A3«<A3+1 . >**0.5 
A3-A3-1 . 

X2iIT)-(l./Al)*A2*A3 
IF(X2 ( IT> . LT. 1 . D-5 )  GOTO  400 
IF( ABS ( CHECK-X2 ( IT ) ) . LE. 1 . D-5 >  GO  TO  400 
300  CONTINUE 

TYPE  * , ' X2  VALUE  HAS  NOT  CONVERGED . ' 
£********★*******************************•*★************************ 
C  WE  NOW  CALCULATE  THE  1ST.  DEGREE  OF  IONISATION  FOR 

C  A  GIVEN  TEMPERATURE. 

£****************************************************************** 
400  Xl(IT)-Al*X2(lT)/2. 

Q*  If*********  ********************************  It*************  ********** 

C  CALCULATE  NATURAL  LOG  OF  THE  TEMPERATURE 

£***************************■**********★***********»*★***********«** 
t(it)-log(sngl(temp) ) 

1000  CONTINUE 

c****************************************************************** 

C  OUTPUT  DATA  TO  A  FILE  FOR  PLOTTING  PURPOSES 

C  OUTPUT  THE  X  AXIS  LABEL 

£****************************************************************** 
WRITE (  2 , 6  0  0 ) 

600  FORMAT!'  NATURAL  LOG  OF  TEMPERATURE’  ) 

C  OUTPUT  A  -l’  FOR  THE  PLOT  FILE  FORMAT 

WRITE! 2 , 700 ) 

700  FORMAT!'  1') 

C  OUTPUT  THE  Y  AXIS  LABEL 

WRITE! 2 , BO  0 ) 
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800  FORMAT (  '  DEGREE  OF  IONISATION  ') 

WRITE(2,?00) 

C  OUTPUT  THE  SUBTITLE 

WRITE! 2 , 960 ) PRESS 
960  FORMAT! ’  PRESSURE-' ,D10 . 3) 

WRITE! 2,950) -NUM 
WRITE! 2 ,900) ATOM 
900  FORMAT!  IX,  6A4) 

C  OUTPUT  NEGATIVE  VALUE  OF  THE  NUMBER  OF  POINTS  TO  BE  PLOTTED 

950  FORMAT! IX, 14 ) 

DO  2000  LK-1 , NUM 

WRITE! 2, 1050 )T(LK) , SNGLtXl (LK) ) 

1050  FORMAT (1X,F6.3,',#,F6.3> 

2000  CONTINUE 

WRITE (2,950) -NUM 
DO  2100  LK=>  1  ,  NUM 
WRITE! 2,1050)t( LK ) ,SNGL<X2(LK) ) 

2100  CONTINUE 

CLOSE! UNIT-1 ) 

CLOSE(UNIT-2 ) 

TYPE  * , *  DO  YOU  WISH  TO  HAVE  ANOTHER  RUN?  ANSWER  Y/N' 

ACCEPT  2110, ANS 
2110  FORMAT! A1 ) 

IF! ANS . EQ . ' N' )  GOTO  3000 
GOTO  1 

3000  STOP  'THATS  ALL  FOLKS' 

END 
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FIGURE  A-l 


Flow  chart  lor  program  used  to  calculate  1st  and  2nd  degree 
ionisation . 
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